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ABSTRACT 

We model the chemical evolution of six UFDs: Bootes I, Canes Venatici II, Coma Berenices, Her¬ 
cules, Leo IV and Ursa Major I, based on their recently determined star formation histories. We show 
that two single-age bursts cannot explain the observed [a/Fe] vs [Fe/H] distribution in these galaxies 
and that some self-enrichment is required within the first burst. An alternative scenario is modelled, 
in which star formation is continuous except for short interruptions when one or more supernovae tem¬ 
porarily blow the dense gas out from the centre of the system. This model allows for self-enrichment 
and can reproduce the chemical abundances of the UFDs in which the second burst is only a trace 
population. We conclude that the most likely star formation history is one or two extended periods of 
star formation, with the first burst lasting for at least 100 Myr. As found in earlier work, the observed 
properties of UFDs can be explained by formation at a low mass (Mvir ~ 10^ Mq), rather than being 
stripped remnants of much larger systems. 

Subject headings: dark ages, reionization, first stars—galaxies: abundances—galaxies: dwarf— 
galaxies: formation—galaxies: star formation 


1. INTRODUCTION 

The discovery of very faint low-mass galaxies, known 
as ultra-faint dwarfs (L < 10^ Lq) has provided a low- 
redshift method of investigating the formation and evo¬ 
lution of the early baryonic systems. UFDs contain only 
old, metal-poor stellar populations and may retain rela¬ 
tively unpolluted chemical signatures of the first gener¬ 
ations of stars. Recent work has ma de progress in de¬ 


termining the star formation history (Weisz et al. 
Brown et al.||2014|) and chemical abundances jFrebel S 


_ II MI _^_ 

Bromm||2(J12||GiImore et al.||2013| [Vargas et al!||2Uld|) of' 

these systems. 

For most stars in UFDs, only medium-resolution spec¬ 
troscopy is available, however this is sufficient to deter¬ 


mine iron and a element abundancesj Kirby et al. | 2009, 
20101. This allows use of the idea of 1 insley (|l979p, who 
suggested that the enhanced [a/Fe] ratio in halo stars 
could be explained by the time delay between Type II 
and Type la supernovae. The progenitors of Type la 
supernovae are evolved low and intermediate-mass stars 
with longer lifetimes than the massive stars that are the 
progenitors of Type II supernovae, such that there is a 
^ 100 Myr period after star formation commences in 
a galaxy when only Type II supernovae enrich the gas. 
Type II supernovae eject much more a elements relative 
to iron than Type la supernovae, resulting in enhanced 
[a/Fe] for stars that form in the first 100 Myr of star 
formation in a system. A galaxy that only contains stars 
with [a/Fe] above or close to the mean of 0.35 expected 
for Type II supernovae in a typical initial mass function 
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(e.g. Salpeter 1955 Kroupa 2001|) is likely to have formed 
stars tor less than 100 Myr ([Gayrel et al.|[2004 Frebel & 


Bromm 


2012 ). 


Enrichment from Type la supernovae is not the only 
possible explanation for low [a/Fe], which can also be 
explained by an initial mass function (IMF) that favours 
lower-mass Type II supernovae, or by types of supernovae 
in which most a elements do not esca pe into the inter- 

Observationally 


stellar medium (Karlsson et al. 
disentangling the contribution o 


2012 ). 
the decline in 


a/Fe] as 

a result of Type II supernovae as compared to Type la 
supernovae is therefore difficult and is likely to require 
high-resolution observations of iron-peak and s-process 
elements. It should also be noted that there is not neces¬ 
sarily a direct relationship between [Fe/H] and time, as 
the merging of gas clouds can create multiple metallicity 
populations ( Wise et ar]|2012 ). 

The low luminosity of the UFDs means that chemical 
abundance data is available for only a few stars in each 
galaxy, resulting in large uncertainties for individ ual sys¬ 
tems. vHowever, the situation is improving, with |Brown| 


et al. (2014) pro viding [Fe/H] for a l arge sample of stars 
m SIX UFDs and Vargas et al. (2013) determining [a/Fe] 


abundances for 61 stars in eight UFDs. The [Fe/H] dis¬ 
tribution from these observations was then used to de¬ 
termine the ages of the stars relative to the M92 globular 
cluster. 


Brown et al. (2014) used isochrone fitting to determine 
the star formation history (SFH) of six UFDs; Bootes 1, 
Goma Berenices, Ganes Venatici H, Hercules, Leo IV and 
Ursa Major 1. Without the constraint of spectroscopic 
abundances for [a/Fe], they found that the SFH could 
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be fit by a two single-age burst model with three param¬ 
eters: the ages of the two components and the proportion 
of stars in each burst. Adding parameters for the dura¬ 
tion of the two bursts did not improve the fit, which the 
authors suggest indicates a narrow age range for the stars 
within each burst. 

In this work we simulate chemical abundances given 
two possible star formation histories us ing the models 
of Mvir = 10^ Mf7) halos presen t ed in Webster et al. 


2.1.1. Bootes I 


(2014) & Bland-Hawthorn et al. (2015). Collins et al. 


(2014) showed that Bootes 1 and Hercules have circu- 
lar velocities of 5-7 kms“^, consistent with halo masses 
this low. We al so investigate whether observed chemi¬ 
cal abundances ([Gilmore et al.||2013 Vargas et al.||2013[ 


Brown et al. 20i4p are consistent with the Brown et al. 
( |2014[ ) star formation histories. 

2. MODELS 

The simulations used to model the chemical abun- 


Hawthorn et al. (20151).^ 

I’he 3U hydro/ionization code 

I'yris Alpha (^Sutherland 

2010 

) was used to model the 


0 

halo. The density and metallicity distribution of the gas 
after the first supernova was then used as a template 
to simulate the effects of later supernovae, allowing an 
estimation of star formation and gas enrichment over pe- 
riods of up to 600 Myr. Using the method of |Argast| 


et al. (2000), a number of cells were randomly selected, 
with stars forming in these cells with a probability pro¬ 
portional to the square of the density of the gas. 

In our model of two single-age bursts, the first burst 
stars form in gas enriched from [Fe/H] = — 4 by only a 
single 25 Mq star. The density and metallicity distribu¬ 
tion for the first burst gas is taken from the hydrody- 
namical model 15 Myr after the supernova. At this time 
most of the enriched gas has returned to the centre of the 
galaxy, such that the number of cells with nn > 10 cm“^ 
is > 70% of the number in the undisturbed state. The 
second burst forms in gas enriched by Type II and Type 


Woosley & Weaver 

(1995 

)& 

Iwamoto et al. 

(1999 


as m 

The number of Type la supernovae is similar to the num¬ 
ber of Type II supernovae. Bootes I is enriched by only 
Type II sup ernovae, because the two bursts in |Brown| 
et al. (2014) are onl y 100 Myr apar t . Th e continuous 
model IS described in Webster et al. (2014). The gas is 
enriched by only Type 11 supernovae tor the first 100 Myr, 
after which Type la supernovae occur with a rate as in 
Jimenez et al. (2014). Because higher mass stars yield 
more alpha elements, varying the mass of the initial star 
would alter [a/Fe] at low [Fe/H], meaning that we should 
not necessarily expect the models to fit the observations 
at low [Fe/H]. However, this effect is washed out after a 
few supernovae. 

2.1. Single-age bursts 

Base d on colour-magnitude diagrams, [Brown et al. 


(2014) modelled the star formation histories of six UFDs 


as two single-age populations, determining the ages of 
the bursts and the proportion of stars in each burst. A 
model with extra parameters for the duration of star for¬ 
mation in each burst did not improve their fit, suggesting 
a narrow age range within each burst. 


The two-burst model of Brown et al. (2014) has 97% 
of the stars in Bootes I forming in the second burst, with 
the two bursts being only 100 Myr apart. This is the only 
galaxy for which the Brown et al. (2014) model favours 
such closely spaced bursts, as well as the only case where 
significantly more stars form in the second burst. 26 of 
the 38 observed stars (68%) in Bootes I have [Fe/H] > 
—2.5, compared to < 50% for the other hve galaxies. Our 
model for Bootes lis shown in the first pa nel of Figs. | T 
and along with the [Fe/ H] data from Br own e t ah 
(201^ and [a/Fe] data from Gilmore et m '[20131). In 
our model, the 3% of stars that form in the hrst burst 
enrich the gas to [Fe/H] ~ —3, with a spread in [a/Fe] 
resulting from variations in the extent to which these 
stars enrich different regions of the galaxy. 

A cross-correlation between the observed [Fe/H] his¬ 
togram and our modelled histogram gives a lag of 0.6 dex 
for the model. This could be explained by a higher initial 
metallicity in Bootes I than in our model, or a top-heavy 
IMF such that a greater proportion of the stars in the 
first burst produced supernovae. The red dotted his¬ 
togram, corresponding to an increased Type H SN rate, 
shows a much smaller lag compared to the obse rvations. 

[a/Fe] observations from Gilmore et al. (2013) indicate 
at least some self-enrichment, with evidence of a decline 
in [a/Fe] with increasing [Fe/H]. This could be explained 
by Type la supernovae, but can also be explained by 
lower mass Type H supernovae, with stochastic sampling 
of a Kroupa IMF selecting many 8-15 Mq star in Bootes 
I, reducing [g/Fe] , If Type H superii ovae are the cause 
of the decline, the Brown et al. (2014) scenario with two 
single-age bursts 100 Myr apart is possible for Bootes I. 

2.1.2. Canes Venatici II 


The Brown et al. (2014) model found that 95% of the 
stars in Ganes Venatici 11 formed in the first burst. Our 
modelled [a/Fe] and [Fe/H] f or this scenario is shown 
along with observational data (Vargas et al.|2013 Brown 
et al. 2014) in the second panel of Figs. T] & Whife 
our model of the two-burst scenario can nt most of the 
stars, it produces insufficient scatter to fit the stars at 
[Fe/H] ^ —1.2. Furthermore, four of the eight stars for 
which [a/Fe] data is available show [a/Fe] suppressed by 
0.5 dex compared to the other four stars. If low [a/Fe] 
reflects self-enrichment, these stars must have formed in 


the second burst. However, Brown et al. (2014) suggest 
that 95% of the stars in Ganes Venatici 11 formed in the 
first burst. Selecting eight stars randomly from such a 
distribution gives a probability of 3 x 10“"^ that four or 
more will be second burst stars. Removing one of the 
high-metallicity stars for which the membership of the 
system is uncertain (J. Simon, private comm.) increases 
this probability to 4 x 10“^. As with Bootes I, the cross¬ 
correlation showed a lag of 0.6 dex in [Fe/H] for the model 
compared to the observations. 

2.1.3. Coma Berenices 

Goma Berenices shows a lag of 0.4 dex in [Fe/H] for 
the best-fittin g m odel compared to the observations. As 
shown in Fig.^ two of the nine stars for which [a/Fe] 
is available show suppressed [a/Fe] abundances, suggest¬ 
ing that they belong to the second burst. Using the same 
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Fig. 1.— f a/Fel vs [Fe/H] for 6 UFDs from our model (black points) and observations (coloured squares) from|Gilmore et al., , . 
(Bootes I) &;|Vargas et al.|(|2013f (the other five galaxies). Dotted error bars represent stars whose membership of the system is m doubt 
(J. Simon, private comm.J 


method as for Canes Venati ci II, this is reasona bly un¬ 
likely {p = 0.05) given that [Brown et al. (2014) deter¬ 
mine that 96% of the stars belong to the first burst. If 
the two-burst model is to explain Coma Berenices and 
Canes Venatici II, it requires more stars in the second 
burst, or enrichment within the first burst. 

2.1.4. Hercules 


insufficient data to conclude whether Leo IV is consistent 


The Brown et al. (|2014 ) model for Hercules has 82% of 
the stars in the hrst burst. As with the previous systems, 
the model overestimates the number of low-metallicity 
stars, with the cross-correlation showing a lag of 0.6 dex 
in [Fe/H] for the model. Fig. shows that [a/Fe] ob¬ 
served in Hercules is lower than predicted by the model, 
suggesting that the gas that formed the first burst of stars 
was enriched by stars with lower alpha abundances than 
a 25 Mq supernova. Five of the 13 stars show sub-solar 
[a/Fe], suggesting Type la enrichme nt. This is possible , 
but unlikely (p = 0.07) given that the Brown et al. (20141 
model has 82% of the stars in the firk burst. 

2.1.5. Leo IV 

Limited abundance data is available for Leo IV, which 
conta ins only four stars wi th known [g/Fe] ( Vargas et al. 
2013) and 13 with [Fe/H] (Brown et al.||2014|). 4'he four 
stars with [a/Fe] abundances are consistent with no de¬ 
cline with increasing [Fe/H] or a slight decline. There is 


with the Brown et al. (2014) two-burst model. 


2.1.6. Ursa Major I 

Ursa Major I has half-light mass of 2.6/)]x 10^ Mg 


(Collins et al. 2014) and a luminosity of 1.4±0.4x 10^ L©, 
which IS much more massive and luminous than our 
model and therefore the results from the model should 
be treated with cau tion. It is the only galaxy studied by 
(2014) with similar numl^rs of stars form- 
3urst. As shown in Fig. the model can 


Brown et al. 
ing in each 


explain the observed abundances except for the stars 
with [Fe/H] ~ —1.5. However , three of these hig her- 
metal licity stars are inclu ded in Vargas et al. (2013) but 
not injBrown et al. (2014) because they are no longer be¬ 
lieved to~5e%nemEers~(Tr Simon, private comm.). When 
these three stars are excluded, the model is a much better 
fit to Ursa Major 1. 

2.1.7. Summary 

The two single-age burst model of|Brown et al. (2014) 
found that four of the six UFDs studied formed 75-967o of 
stars in the first burst. However, 38% (13/34) of stars in 
these systems show [a/Fe] suppressed >0.5 dex compared 
to the high [a/Fe], low [Fe/H] stars in the same galaxy. 
These stars show [a/Fe] <0.2, indicative of enrichment 
either from Type la supernovae or lower mass Type H 
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[Fe/H] [Fe/H] [Fe/H] 


Fig. 2.— Histograms of [Fe/H] from our model (green) and the iBrown et al.| (|2014j observations (blue) for 6 UFDs. The modelled 
histogram is normalised to the number of observed stars in each Uh L). U he dotted histograms correspond to a model with a higher Type 
II SN rate and lower Type la SN rate (red), and a lower Type H SN rate and higher Type la SN rate (light blue). 


supernovae after the time at which the high [ct/Fe] stars 
formed. To fit observations of [ct/Fe], the two single-age 
bursts model requires a more even distribution of stars 
between the first and the second burst. 


2.2. Continuous model 


In Webster et al. (2014|, we presented a star formation 
history where each supernova in a 10^ Mq galaxy tem¬ 
porarily blows out the gas from the centre of the system, 
pausing star formation for ~ 10 — 15 Myr until the gas 
recovers. Star formation proceeds in the intervals be¬ 
tween the supernovae, which are typically ~10-20 Myr, 
but can be as long as 30-40 Myr. Longer, possibly per¬ 
manent pauses may be caused by a large number of su¬ 
pernovae occurring nearly simultaneously. The output 
from this model is present ed in Figs. fc [4| al o ng with 
observed ab u ndances from [Brown e t m. (2014); Vargas 
et al. (2013); Gilmore et al.] 1 2Dl3p . 


Like the two-burst model, the continuous model shows 
a clump of stars at low [Fe/H] with [a/Fe] ^0.8, how¬ 
ever this is a smaller proportion of stars than for the 
bursts model. The gas is then enriched by Type II su¬ 
pernovae, resulting in stars forming with [Fe/H] ^ — 2.5 
and [a/Fe] ~ 0.6. This level of [a/Fe] is higher than the 
average of ~ 0.35 from Type H supernovae because the 
low star formation rate results in a low supernova rate, 
meaning that the enhanced a abundances present in the 


initial gas and from the 25 Mq star have not yet been 
washed out. After 100 Myr, Type la supernovae enrich 
the gas, eventually reducing [a/Fe] to sub-solar levels. 

The [Fe/H] histogram for this model produces signifi¬ 
cantly more higher metallicity stars than the two single¬ 
age bursts due to the self-enrichment resulting from ex¬ 
tended star formation. The number of stars at each 
metallicity increases until [Fe/H] = —2.5, as expected for 
a reasonably constant rate of enrichment, because [Fe/H] 
is a logarithmic scale. There is then a decline caused by 
the start of Type la enrichment. Type la supernovae 
yield much more Fe than Type H, resulting in a more 
rapid increase in [Fe/H] and therefore fewer stars at a 
given metallicity. The number of stars at each [Fe/H] 
then begins to increase again, followed by a decline due 
to the truncation of star formation. While the observed 
[Fe/H] histograms contain limited stars, there may be 
signs of this pattern in Coma Berenices, Ursa Major I & 
Hercules, all of which show valleys near [Fe/H] = — 2.5, 
which is the metallicity at which these systems start to 
show evidence of Type la supernovae. 


2.3. Comparison between the models 

To compare the models we implemen t a me thod based 
on Section 3.2 of Price-Whelan et al. (2014|). For each 
point in the model, the probability of a modelled star k 
matching the observed star s is: 
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Fig. 3.- 
squares). Dotted 


/Fe] vs [Fe/H] in our continuous model (black points) and as observed (|Gilmore et al.j|2013| |Vargas et al.||2013|| (coloured 
:ed error bars represent stars whose membership of the system is in doubt (J. bimon, private comm.J. 


Pks — 


2,7 T CxC y 


/ / (^mod ^obs) , (ymod Yobs) \\ 

exp(-(-- ■)) 




2a2 


d) 

where x = [Fe/H] and y = [a/Fej. This is then summed 
over all the modelled stars and normalised by the num¬ 
ber of stars in the model. The overall likelihood for the 
system is then the product of the likelihoods of the indi¬ 
vidual observed stars: 


L = \ls{^Y.pks) ( 2 ) 

This gives the log likelihoods shown in Table 1. In 
nearly all cases the extended model has a greater likeli¬ 
hood than the two-burst model, with the only exception 
being Ursa Major I when the stars with questionable 
membership are excluded. 

3. CONCLUSIONS 

We have modelled the chemical evolution resulting 
from two possible star formation histories. Our conclu¬ 
sions are as follows: 


1. The t wo single-age burst model of Brown et al.| 
(20141 produces too many stars with low [Fe/HJ, 
with a lag of ~ 0.6 dex in the model compared 
to the observations. This could be explained in 


TABLE 1 

Log likelihoods of the models. The numbers in brackets 

ARE THOSE OBTAINED WHEN STARS THAT ARE NOW BELIEVED NOT 
TO BE MEMBERS OE THE SYSTEMS ARE EXCLUDED (J. SiMON, 
PRIVATE COMM.). 


System 

Two-burst 

Extended 

Bootes I 

-12.0 

-6.8 

Canes Venatici II 

-20.3 (-14.4) 

-4.2 (-3.2) 

Coma Berenices 

-4.9 

-3.6 

Hercules 

-14.2 

-9.4 

Leo IV 

-3.0 

-2.3 

Ursa Major I 

-15.3 (-3.5) 

-8.2 (-4.2) 


part by the systems forming at a higher metallicity 
[Fe/H]~ —3.5, rather than —4 as assumed in our 
model. 


2. Extended star formation is a better fit to obser¬ 
vation s of [a/Fe] in UFDs than the Brown et al.j 
(20141 two single-age bursts model tor all galax- 
ies except Ursa Major I. The difference is largest 
for Canes Venatici II, in which half the observed 


stars show suppressed [a/Fe], while the Brown 


et al. (20141 star formation history has 95% ot stars 
forming m the first burst. 


3. Enrichment within a burst is required to explain 
the number of low [a/Fe] stars in systems where 
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[Fe/H] [Fe/H] [Fe/H] 


Fig. 4.— [Fe/H] histograms for stars in our continuous model (blue) and the [Brown et al.| ||2014|l observations (green). The dotted 
histograms correspond to a model with a higher Type II SN rate and lower Type la SiN rate (red), and a lower Type II SN rate and higher 
Type la SN rate (light blue). 


Brown et al. (2014) determine that the vast ma- 
jority of stars are in the first burst. The suc¬ 
cess of the Ursa Major I two-burst model suggests 
that two-burst models can explain at least some 
UFDs, however a more even distribution of stars 
between the two bursts is required. An alternate 
explanation is extended bursts, such that the first 
stars in each burst enrich those formed later, ei¬ 
ther through Type la supernova enrichment, with 
a timescale >100 Myr, or lower mass Type II su¬ 
pernovae (> 25 Myr). 

4. Our modelled Mvir = 10^ Mq systems with ex¬ 
tended star formation can reproduce the observed 
chemical abu ndances. This provide s support to the 
conclusion of Webster et al. (2014) that the UFDs 
can be explained as systems with a low forma¬ 
tion mass rather than stripping from much larger 
(Mtot 10® Mq) halos. 


determination of the abundances of more elements. This 


A co mbination of isochrone fitting as in Brown et al. 
(2014), observations of chemical abundances, and mod¬ 
elling of chemical evolution can reveal the star formation 
history of UFDs. Knowledge of the star formation histo¬ 
ries can give insight into star formation processes prior 
to the epoch of reionization. Deeper spectroscopy with 
the next generation of extremely large telescopes will al¬ 
low the study of a larger sample of stars, as well as the 


should allow chemical tag ging (Bland-Hawthorn et al. 
20101 iKarlsson et al.||2012 ). If stars in dwarf spheroidals 
and Uh'Ds formed in large clusters, it should be possible 
to identify stars born within the same cluster, as there 
will b e very low scatter in c hemical abundances between 
them. Karlsson et al. (2012|) presented tentative evidence 
for a cluster in Sextans at [Fe/H] = —2.1 based on three 
stars with similar [Fe/H] abundances that also had simi¬ 
lar Mg, Ti, Cr & Ba abundances. An equivalent study of 
UFDs to identify individual clusters requires more chem- 
ic al abundance observat ions than is currently available. 


Karlsson et al. (2012) also found that the cumulative 
metallicity functions of UFDs showed less clustering than 
for dSphs. They suggest that this could result from 
UFDs forming before reionization at masses bel ow the 
atomic hydrogen cooling limit M ~ 10® Mq as in Bovill 
& Ricotti (2009). Star formation in such halos would be 
affected by inefficient cooling and feedback from Lyman- 
Werner radiation, which dissociates H 2 . This could result 
in a lower cluster mass. Lower cluster masses provide 
support for the result in this paper that UFDs experi¬ 
enced more extended star formation than in the case of 
large clusters with multiple supernovae. 
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